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Abstract 

We represent the free energy functional by a diagrammatic series with tensorial coefficients 
indexed by powers of length scale. For hard cores, we obtain Percus' exact functional in one 
dimension and the Kierlik-Rosinberg form of fundamental measures theory in three dimensions. 
In five dimensions, the functional describes bulk fluids better than Percus- Yevick theory does. 
At planar walls density profiles oscillate with smaller periods than in lower dimensions. Our 
findings open up avenues for treating both more general high-dimensional systems, as well as 
three-dimensional mixtures via dimensional reduction. 
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Studying the effects of inter-particle hard core repulsion has provided significant insights 



into interfacial and capillary (phase) behaviour of simple and complex liquids Recent 
examples include studies of ion-specific excluded- volume correlations and solvation forces 2|, 
and the wetting properties of a solid substrate by a "civilized" model of ionic solutions 3|. 
Density functional theory (DFT) {4] is a primary tool for the investigation of such situations. 
Applying DFT requires one to have an approximation for the central (and in general un- 
known) Helmholtz free energy as a functional of the one-body density distributions Pi{r) of 
all species r is the spatial coordinate. One such example, which has become a cornerstone 
of DFT, is Rosenfeld's fundamental-measure theory (FMT) 5|-[8| for additive hard sphere 
mixtures. 

There is considerable current interest in the study of the hard sphere model in large 
space dimensions D. Motivation for such work originates both from a desire to develop and 
to gain new insights into the general structure of liquid state theories, as well as from the 
question whether anything is special about D = 3. Studies that exemplify this strategy 



include investigations of fluid structure |9l-ll2| and of freezing 13 



Much less is known 



about inhomogeneous hypersphere fluids. In this Letter we present a formalized framework 
for a generic hard hypersphere density functional. Its structure is based on a density series, 
where the central approximation lies in the star-like topology of the integral kernels that 
couple the field points together [cf. Eq.([3])]. 

The lowest order contribution to the exact virial expansion of the excess (over ideal 
gas) free energy functional is —ksT f dripi{ri) f dr2Pj{r2) fij{\ri — r2|)/2, where the Mayer 
bond fij{r) for hard spheres as a function of distance r is fij{r) = —1, if the two spheres 
with radii Ri and Rj overlap, and zero otherwise; ks is the Boltzmann constant and T is 
temperature. Graphically this can be represented by two filled dots (each one corresponding 
to one multiplication by density and one integration), that are joined by a line which stands 
for fij{r). The Kierlik-Rosinberg (KR) deconvolution of the Mayer bond for arbitrary (odd) 
D is /jj(|ri— r2|) = — J2^=o I '^^'^lli\'''^'^~^\)^D-^li\'''^'i~^\)^ with species-dependent geometric 
"weight" functions w^(r); all Greek indices run from to D here and in the following. Hence 
diagrammatically 

ij ij fi,u=0 
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where the cross indicates the position x with integration and the bonds labelled /i and 
V represent ti'^dri — x|) and w^(|r2 — x|), respectively. The second-rank tensor M has 



components M^^ = M^^, = 1 for fi + u = D and is zero otherwise [IGJ. The third diagram 
in ([T]) is a short-hand notation, where e.g. the first bond, together with its filled circle, 
represents Xlj / dripi{ri)w^^{\ri — x|) = ri^(x); tensor contraction over all Greek indices is 
implied. Both and n,, possess units of (length) In D = 3, the n^(x) are the familiar 
KR weighted densities p, if the corresponding form of the w^^Jr) is used. 

We approximate the exact "triangle diagram" of third order 17|] by a three-arm star: 
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where the third-rank tensor J has components J^'^^, which we assume i) to be symmetric 
under exchange of indices, ii) to satisfy J"^'"^ = M^'^ , and iii) to vanish ii ii + v + t ^ 2D, 
because the diagram has to be of dimension (length) °. 

Restricting ourselves to star topology, we assume the excess free energy functional in D 
dimensions to have the generic form 



/3J^exc[{Pi}] = - 



M 

-X- 




+ 



(3) 



where /3 = 1/(/cbT), and the scalar coefficients are taken from the zero-dimensional excess 
free energy V5o(^) = (1 ~ "'?) —''?)+ ^ = Yl'k=2 /[k{k — 1)] = rf (2 + rf /Q + rj^ /12 H — • ; 
here ?7 is a dummy variable. The product JJ indicates tensor contraction via X]!^=o J^^^ Jt'^^^ 
where J^"^^ = Y1!!!'=q^tt'J^'^^'-, the triplet JJJ represents X]fo-=o J J J ^tc. In order 
to find a closed expression for ([3]), we define a matrix of weighted densities 



N(x)= |. 



(4) 



where the spatial argument x is not integrated over. The components of N(x) are N^^{x) = 
'^i'^^=Qjf^^^{wl * Pi)(x), where * denotes the D-dimensional convolution {g * /i)(x) = 
J drg{-x — r)h{r) of two functions g^x) and /i(x). The definition (jl]) allows one to express 



the k-th order diagram in ([3]) as a matrix power N(x) with i nteg ration over x. Subsequently 



one has to take the appropriate component n = 0,1/ = D [l8|. Hence we can view ([3]) as 
one component of the tensorial functional 

^^E^(^ = J d^M^i^)), (5) 

where the second equality comes from using the Taylor series of ^poiv)- The physically 
relevant functional J-'cxc[{pi}] is the ji = = D component, i.e. J-'exc = (Fcxc)o^) which 
can be written as /3J-'exc[{pi}] = / rfx$, where the free energy density $ = [(y9o(N(x))]Q'^. 

All properties specific to the value of D are entirely encapsulated i) in the dimensionality 
of the space integrals, ii) in the prescription of the D + 1 weight functions if^(r), and iii) in 
the choice of the (constant) components J^'^'^ of J. For all D > 0, we use ^^^(r) = 0(r — Ri) 
and w^Q_i{r) = 6{r — Ri), where 9(-) is the Heaviside (unit step) function and 6{-) is the 
Dirac delta distribution. For all cases considered, we could obtain the specific form of the 
further weight functions and the form of J from the minimalistic requirement that the two- 
body bulk fluid direct correlation function Cy(|r — r'|) = — 5^/3J^p^j,/5pj(r)5pj(r')|p^,=const is 
free of divergences, i.e., 

|cij(^)| < oo, for all r. (6) 
In D = 1, J°^^ = 1 determines J entirely and straightforward algebra yields PTcxc = 



J (ixr?,o(x) ln(l — ni{x)), which is exact 19|, |20|. In D = 3 we use the KR form for the 



additional weight functions w\{r) and WQ{r), and obtain J^^^ = J^^^ = 1, J^^^ = l/(47r), 
which gives the KR functional [6] . In D = 5 we obtain 

nin4 + n2n3 n2nl + nln^/ {Gi-n'^) nsnl 
$ = -no ln(l - ng) H \ h ,^ r^r + 



l-ris (1-^5)2 487r2(l-n5)3 160712(1-^5)4 

(7) 

with the weight functions = Q{Ri — r), wl = 6{Ri — r), = 6', w\ = 6' / (IGir'^r), 
w{ = (8r-i(5"-5"')/(647r2), = {2Ar~^6' + 2Ar-^6" -12r-^6"' + 6"")/{6A7r^), where 6'{Ri- 
r), . . . , 6""{Ri — r) denote successive derivatives of the Dirac delta function; the arguments 
have been left away for clarity. 

For constant density fields the 5D functional yields an analytical excess free energy for 
bulk fluids, -Fexc({pi}) = -^cxc[{pi = const}]. When scaled as /3Fexc/V this equals (^^, but 
with = "^iCfiPi = const, where the fundamental measures, = ^ '^rw^(r), are = 1, 



= i?., = = 327r2i?f/3, = 8n^Rf/3, and C5 = Sn'^Rf/lb. The total packing 

fraction is simply 7] = n^. The pressure is obtained from P = —dFc^d dV + ksTp, where 
p = 'Y^j^Pi- Fig. [T] shows PP/p for pure 5D hard hyperspheres as a function of rj. The 
DFT result lies between the two Percus-Yevick (PY) equations of state; these deviate more 
strongly from each other than in three dimensions (3D). Moreover, the DFT only slightly 
underestimates, over the full range of (fluid) packing fractions, the semi-empirical Luban- 
Michels (LM) equation of state [9], which yields excellent agreement with simulation data. 
For two representative binary mixtures (inset of Fig. [1]), the DFT slightly underestimates 



10| , with similar (small) deviations as in the pure 



the very accurate equation of state of Ref. 
case above. 

Fig. [2] shows the pair direct correlation function c(r) for the pure 5D hard hypersphere 
fluid. Outside of the core, the direct correlation function vanishes, c(r > o") = 0, and it is 
a fifth-order polynomials in r for r < a, similar to the PY result; a = 2R is the particle 
diameter. However, the DFT result is smaller in magnitude than in PY theory. This 
constitutes an improvement, as can be observed by comparing to the sophisticated theory 
by Rohrmann and Santos (RS) pjj], which is built on the LM equation of state. For 77 = 0.2, 



a value higher than that at freezing 



13|, we find that indeed the DFT result deviates much 



less from the RS theory than PY does. The inset of Fig. [2] shows an example of partial 
pair direct correlation functions, calculated from the functional. As expected, these satisfy 
Cij{r > Ri+Rj) = and they are free of divergences. Fig. [3] displays corresponding results for 
the radial distribution function g{r), obtained from the Ornstein-Zernike (OZ) relation with 
c(r) as input. Outside of the core the deviation between PY and DFT is very small, but the 
contact value from DFT is lower and some core violation is apparent. We also show results 
from a test particle calculation, where we have numerically minimized the density functional 
in the presence of an external potential that equals the pair interaction potential. Here the 
core condition is automatically satisfied. The contact value g{cr^) increases, and becomes 
larger than the PY value. When input into the virial theorem, PP/ p = 1 + 2^^^r]g{a^), the 
result is almost indistiguishable, see Fig. [H from the quasi-exact LM result. The inset of 
Fig. [3] shows the partial pair distribution functions gij{r) for the mixture considered above 
(inset of Fig. [2]). The agreement with simulation results [12:] is excellent. Hence the DFT 
gives very good account of both the thermodynamics and the structure of pure and binary 
bulk liquids of 5D hard hyperspheres. 
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Density profiles at a hard wall are shown in Fig. Hl^a) for the cases ID, 3D and 5D, 
keeping the scaled bulk density pa^ = 0.7 = const. With increasing dimensionality, the 
distance between the liquid layers that are induced by the wall, decreases strongly. The 
surface tension 7 = (f2 + pV)/A, where Q is the grand potential, A is the area of the wall, 
and V is the system volume 2; > 0, is shown in the inset of Fig. SJ^a) 2l|]. This is seen to rise 
much more strongly with rj for increasing D. Fig. Hl^b) shows density profiles in a slit for 
ID (hard rods on a line), 3D, and for 5D. Again very different structuring is apparent, with 
the position of the secondary peaks being closest to the wall for 5D and five (four) layers 
present in 5D (3D). The inset of Fig. ID^b) shows the solvation force /3/s(L), which for hard 
walls can be obtained from /3/g(L) = p(cr+/2) — /3P, as a function of the slit width L. This 
again highlights sensitive dependence on dimensionality. 

Our formalized diagrammatic approach to DFT has significant relevance for future de- 
velopments: i) A whole class of FMT functionals for 3D models can be generalized in 
a straightforward way to high dimensions. This includes the penetrable sphere model of 
particles interacting with a repulsive step function pair potential [24] (via adapting the 
scalar coefficients in Eq. (|3]) in order to allow for thermally induced overlap of particles), 
the Asakura-Oosawa colloid-polymer mixture j25| (via linearizing the hypersphere mixtures 
functional in one of the densities, which then describes an ideal polymer species), or the 



Widom-Rowlinson model 



26| (by using the appropriate expression for the zero- dimensional 



free energy in order to generate the scalar coefficient in Eq. ([3])). ii) Novel functionals for 3D 
models can be obtained via dimensional reduction, where the hard hypersphere functional 
(or one of its descendants above) is exposed to dimensional crossover, such that the den- 
sity distribution is assumed to be a delta function (or superposition several thereof) in one 
or more space dimensions. The particles in each hyperspace then constitute an individual 
species in 3D. Mixtures with negative non-additivity can be obtained via this procedure, 
see e.g. Ref. js?] for a corresponding example of crossover from 3D to 2D. iii) Although we 



have used the KR approach, very interesting and useful 



from generalizing Rosenfeld's more geometric formulation 5[ to high dimensions. This also 



'urther insights could be gained 



22 



23|, possibly 



applies to its refinements based on considering zero-dimensional cavities 
along the initial steps taken in Ref. 15 |. 
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FIG. 1: Compressibility factor PP/ p as a function of packing fraction r] for 5D hard hyperspheres, 
as obtained from the DFT bulk free energy (thick line), and compared to the PY compressibility 
(PYc) and virial (PYv) results (dashed lines), and the LM (thin line) equation of state [9(]. Also 
shown are DFT results from the test particle limit (TPL, symbols) using the virial theorem. The 
inset shows PP/ p as a function of the total packing fraction r] = for binary mixtures of 5D hard 
hyperspheres as obtained from the DFT (solid lines), along with the equation of state of Ref. [lol | 
(dashed lines). The upper pair of curves is for size ratio 4 and composition 0.75 of the larger 
species; the lower pair of curves is for an equimolar mixture with size ratio 2.5. 
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FIG. 2: Bulk two-body direct correlation function, c(r), as a function of the scaled distance r/a, 
for 5D hard hyperspheres, as obtained from the DFT via second functional derivative (solid lines), 
along with the PY result (dashed lines), for packing fractions r] = 0.05, 0.1, 0.15, 0.2_(from top to 



bottom). Also shown is the result of the theory by Rohrmann and Santos (RS) [ll[ (symbols), 
which is based on the LM equation of state [9], for r] = 0.2 (we omit the small, positive Yukawa tail 
outside of the core). The inset shows DFT results for the partial pair direct correlation functions 
Cjj(r) as a function of r/ai, where 0"i = 2Ri, for a binary mixture of 5D hard hyperspheres with 
size ratio 4, composition of 0.75 for the (larger) species 2 and total density paf = 1.4, corresponding 
to ?? = 0.173 (cf. inset of Fig. d]). 
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FIG. 3: The radial distribution function, g{r), as a function of r/cr, for 5D hard hyperspheres as 
obtained from the DFT using the OZ equation (sohd hne), using the test particle method (symbols), 
and compared to the PY result (dashed line) for packing fraction r/ = 0.2. The OZ result reaches 
g[^) = —2.82 at zero separation. The inset shows the partial pair distribution functions gij{r) as 
obtained from the OZ route (lines) for the binary mixture described in the inset of Fig. [2l along 
with simulation results jl^ (symbols). 
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FIG. 4: Scaled density profiles p{z)a^ as a function of the scaled distance z/a for packing fraction 
r? = 0.7 in ID (dotted line), 0.367 in 3D (dashed line) and 0.115 in 5D (solid line), (a) Hard 
hypersphere fluids adsorbed at a hard wall; z is the distance from the wall. The inset shows the 
surface tension Pja^~^ due to the presence of the wall as a function of 77. (b) Density distributions 
for hard hypersphere fluids in a planar slit with width L = 4.5(T in chemical equilibrium with a 
bulk having the same value of r] as given above. The inset shows the scaled solvation force Pfs^^ 
as a function oi L/a. 
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